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Resumen 


La estimación de los requerimientos de agua de los cultivos usando el 
método dual de FAO-56 necesita de la parametrización del coeficiente 


del cultivo (XK) en forma robusta y confiable. El uso de sensores 
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remotos, a través de índices de vegetación (IV), para definir una 
relación con K¿ es una tarea que debe ser entendida en función de sus 
alcances y limitaciones. En este trabajo se pone en perspectiva qué es 
lo que estiman los IV en relación con los Kc, definiendo que su uso está 
limitado a no considerar el efecto de la humedad del suelo en estos 
últimos. Usando el índice IV_CIMAS y su aplicación a dos experimentos 
(maíz y trigo) con mediciones de evapotranspiración en lisímetros se 
analizan las relaciones de índice con Kc, a fin de establecer una 
metodología para parametrizar sus relaciones en forma sólida, usando 
los inicios y finales de las etapas de desarrollo de los cultivos por medio 
de un modelo expo-lineal del crecimiento. Las relaciones obtenidas en 


la parametrización resultaron en buenos ajustes (R? > 0.96). 


Palabras clave: KIMO SAVI, medio radiativo equivalente, maíz y 
trigo, coeficiente de cultivo. 


Abstract 


The estimation of water requirements of crops using the FAO-56 dual 
method needs a robust and reliable parameterization of the crop 
coefficient (Kc). The use of remote sensing through vegetation indices 
(VI) to define a relationship with Kc is a task that requires an 
understanding of the scope and limitations of each. This paper sets the 
perspective about estimates based on VI relative to those based on Kc, 
noting that the latteris rather limited for not considering soil water 
content. Using the KIMO SAVI index and its application to two 
experiments  —(maize and  wheat) with  evapotranspiration 
measurements using lysimeters, the relationship of this index with K. 


is analyzed aiming to establish a methodology for a solid 
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parametrization of the relationships between these two approaches, 
using the beginning and end of the crop development stages with an 
expo-linear growth model. The relationships obtained in this 
parameterization yielded a good fit (R? >0.96). 


Keywords: KIMO SAVI, equivalent radiative media, com and wheat, 
crop coefficient. 


Recibido: 05/08/2020 
Aceptado: 04/05/2021 


Introducción 


La asistencia satelital en nego usando productos amigables es una de 
las tareas actuales de mayor interés en las aplicaciones de la tecnología 
de los sensores remotos (Calera, Jochum, Cuesta, Montoro, € Fuster, 
2005). La problemática asociada con las estimaciones satelitales de la 
evapotranspiración actual de los cultivos (ET) para definir los 
requerimientos de agua de los cultivos necesita enfoques pragmáticos, 


pero sólidamente fundamentados en teoría y evidencia experimental. 


La tecnología de los sensores remotos permite el seguimiento del 
crecimiento y desarrollo de la vegetación a través de las relaciones 


teóricas o empíricas estimadas entre las variables biofísicas de los 
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cultivos y los índices espectrales de vegetación (IV), de los cuales 
existe una gran cantidad (Verstraete €: Pinty, 1996; Paz et al., 2015). 

La gran mayoría de los IV están basados en el alto contraste 
entre la banda del rojo (R) y la del infrarrojo cercano (IRC) para la 
vegetación (Tucker, 1979). Los índices de vegetación están 
influenciados por múltiples factores, incluyendo los efectos 
atmosféricos, la geometría sol-sensory la mezcla suelo-vegetación. En 
particular, el efecto del suelo (fondo de la vegetación) es relevante 
para estimar las variables biofísicas y el coeficiente del cultivo. Se han 
utilizado diferentes IV que minimizan el efecto del suelo, tales como 
los índices GESAVI (Gilabert, González-Piqueras, García-Haro, 8 Meliá, 
2002) y NDVIcp (Paz et al., 2007), cuyo rango de validez está limitado 
a la primera fase exponencial de los patrones espectrales asociados 
con el crecimiento de la vegetación (líneas iso-IAF). Romero et al. 
(2009) generalizaron el NDVI¿p para el caso de las fases exponencial y 
lineal (hasta el IAF máximo) usando el índice IV_CIMAS. 


La evapotranspiración para un cultivo dado (ET«) se calcula 
multiplicando la evapotranspiración de referencia (pasto o alfalfa) por 
un coeficiente del cultivo (Kc = ETc«/ ETo). La evapotranspiración de 
referencia fue estandarizada en relación con un cultivo virtual con 
características biofísicas y aerodinámicas definidas (Allen, Pereira, 
Raes, € Smith, 1998). 


Partiendo de la suposición de que las variaciones temporales de 
los IV reflejan los patrones temporales de los coeficientes del cultivo 
(Jackson, Pinter, Reginato, € Idso,1980), se ha propuesto el uso de los 
IV para estimar el Kc (Heilman, Heilman, 8, Moore, 1982). Así, se han 


propuesto relaciones lineales entre el K¿ y diferentes IV (Bausch 8 
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Neale, 1987; Bausch € Neale, 1989; Bausch, 1993; Bausch 1995), 
para usarse en la programación de riegos. 


En un contexto general y desde el punto de vista de eficiencia 
hídrica, caracterizar el intercambio de materia y energía entre la 
atmósfera y la vegetación requiere estimar el índice de área foliar 
(TAP), altura (h) y fracción de cobertura de la vegetación (fv). Esto es 
cierto tanto para el uso de modelos de transferencia de dos o más 
fuentes (Shuttleworth £ Wallace, 1985; Choudhury 8 Monteith, 1988) 
como para aproximaciones operacionales del balance de energía como 
el modelo SEBAL o Surface Energy Balance Algorithm for Land 
(Bastiaanssen, Menenti, Feddes, € Holtslag, 1998) y el METRIC o 
Mapping EvapoTranspiration at high Resolution and with Internatized 
Calibration (Allen, Tasumi, € Trezza, 2005), incluyendo los métodos 
simplificados propuestos por FAO (Allen et a/., 1998). 


En este trabajo se analizan los alcances y limitaciones de la 
metodología de usar IV para estimar coeficientes del cultivo, y se 
introduce un esquema operacional para el uso de los IV optimizados 


en las estimaciones de Kc. 


Teoría 
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Índices de vegetación espectrales, variables 
biofísicas y K. 


El análisis de los patrones asociados con los índices espectrales de la 
vegetación (IV) garantiza un marco teórico de referencia que permite 
revisar su congruencia en relación con la teoría de la transferencia 
radiativa en la vegetación. Este formato de análisis evita el uso de 


consideraciones empíricas para definir que IV es mejor. 


En la Figura 1a se muestra el patrón temporal de crecimiento de 
un cultivo en el espacio espectral R-IRC, simbolizado por isolíneas del 
índice de área foliar (iso-IAF), que representan la misma cantidad de 
vegetación para suelos ópticamente diferentes, la cual se generó 
usando seistipos de suelos (S2, S5, S7, S9, S11 y S12; del más oscuro 
al más claro). Paz, Palacios, Mejía, Martínez y Palacios (2005) detallan 
las simulaciones radiativas mostradas en la Figura la. En la Figura 1b 
se muestran los resultados de un experimento de maíz donde se usaron 
charolas con suelo deslizables debajo del cultivo para simular 
propiedades ópticas diferentes de suelos (Bausch, 1993). 
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Figura 1. Patrones espectrales en el espacio del R-IRC. (a) 
Simulaciones radiativas y (b) experimento de maíz. 
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En la Figura 1 se ajustaron líneas rectas (interacciones radiativas 
de primer orden) a los valores iso-IAF de la reflectancia de la 
vegetación, donde las reflectancias del suelo varían (línea del suelo en 
la Figura 1). El análisis de la Figura 1 define varios patrones 
importantes para entender el comportamiento de la reflectancia y de 
los IV durante el desarrollo de los cultivos: a) si se unen los valores de 
igual índice de área foliar (iso-IAF) correspondientes a cada curva de 
igual suelo (iso-suelo), se obtiene un patrón cuasi-lineal. Así, el estado 
de crecimiento de un cultivo (IAF) se refleja en una línea recta, 
independientemente del tipo de suelo de fondo en el cultivo; b) la 
pendiente (bo) e intersección (ao) de las isolíneas del IAF varían con el 
valor del JAF. La pendiente de las isolíneas de JAF parte desde una 
pendiente igual a la de la línea del suelo (IAF = 0) y aumenta en el 
sentido contrario a las manecillas del reloj hasta alcanzar un ángulo de 
90%. Esta última condición corresponde al caso de saturación de la 
reflectancia de la banda del R, representada en la Figura 1 como los 
valores de reflectancia arriba del ápice del “sombrero de tres picos” 
(IAF > 5); c) todas las isolíneas del suelo convergen al mismo punto 
de saturación de las bandas visibles. En realidad, la convergencia es 
hacia una línea recta (ápice en la Figura 1), pues cuando la banda del 
R se satura, el IRC no lo hace y sigue creciendo hasta su propio punto 
de saturación. El punto de saturación, llamado también de reflectanda 
infinita (Roo o JRCoo) o de medio denso, es función de las propiedades 
ópticas de las hojas y de su distribución angular (Ross, 1981). 


En la Figura 2a, caso de las simulaciones radiativas, se muestra 
la relación entre los parámetros ao y bo de las líneas rectas iso-IAF para 
el ciclo de crecimiento del cultivo que va desde el suelo desnudo hasta 


la condición de cubrimiento completo del suelo por la vegetación. Se 
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observa que en la etapa inicial hay un patrón de comportamiento tipo 
exponencial hasta un punto donde la banda R se satura 
(aparentemente). Después del punto de saturación aparente de la 
banda R, el patrón es de tipo lineal. El punto inicial de la curva ao-bo 
representa el caso de suelo desnudo (ao = as y bo = bs). El punto donde 
ao alcanza su valor máximo (punto de transición del patrón exponencial 
al lineal) representa el término de la fase de crecimiento exponencial 
e inicio de la fase lineal de la etapa de crecimiento vegetativo de la 
vegetación (Paz, Casiano, Zarco, €: Bolaños, 2013). El punto final del 
patrón lineal de la curva ao-bo representa la situación donde la banda 
del IRC se puede o no saturar, lo cual ocurre cuando el JAF alcanza su 


valor máximo. 
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Figura 2. Patrones de los parámetros ao y bo de las líneas iso-IAF: a) 
simulaciones radiativas y b) experimento de maíz. 


Tecnología y ciencias del agua, ISSN 2007-2422, 13(3), 374-424. DOI: 10.24850/j-tyca-2022-03-08 


2022, Instituto Mexicano de Tecnología del Agua 
Tecnología y Open Access bajo la licencia CC BY-NC-SA 4.0 


Ciencias¿Agua (https://creativecommons.org/licenses/by-nc- 
sa/4.0/) 


En la Figura 2b, datos experimentales del maíz, el patrón entre 
ao y bo es similar al caso de las simulaciones radiativas, pero con 
diferencias importantes en función de las fases del crecimiento de los 


cultivos. 


En el caso de las simulaciones radiativas, la aproximación de 
medio turbio (Ross, 1981) supone un cubrimiento homogéneo y 
completo del suelo por la vegetación, además de que se usó la hipótesis 
de que las propiedades ópticas y geométricas de las hojas no cambian 
durante el crecimiento (incremento en el 7AF) del cultivo. Bajo estas 
consideraciones, el patrón temporal del TJAF puede modelarse con una 
función expo-lineal que suaviza la transición de la fase exponencial a 
la lineal (Goudriaan € Monteith, 1990) (Figura 3a). La Figura 2a 
muestra que en este caso la fase expo-lineal tiene un patrón similar en 
el espacio ao-bo. El modelo expo-lineal, para la variable V en función 


del tiempo t, está definido como: 


V(0) =7In(1 + explr(t— t,)1) (1) 


Con las siguientes características: 


Cc 0.54-1 
Vlt=tp)=2;5tp=t,+ (2) 


Tr 
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Figura 3. Patrones temporales del IAF: a) simulaciones radiativas y 


b) experimento de maíz. 
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En el modelo expo-lineal (Figura 3a), el punto de transición entre 
la fase exponencial y la lineal o 74F7 tiene un tiempo asociado tr. La 
fase lineal extendida hasta cortar el eje del tiempo (Figura 3a) está 
definida por el "tiempo perdido” o tp. Este tiempo es el que se pierde 
por la etapa exponencial hasta que se alcanza la tasa máxima de 
crecimiento, la cual se mantiene durante la fase lineal. El tiempo tmx 
está asociado con el IAF máximo, JAFmx, de la fase lineal (no es el IAF 
máximo de toda la etapa vegetativa). 


El modelo expo-lineal tiene bases biofísicas para el IAF 
(aproximación de medio homogéneo en términos físicos y de medio 
turbio en términos radiativos) (Ferrandino, 1989). Este modelo ha sido 
aplicado exitosamente en el análisis de la variación temporal de los 
índices de vegetación (Odi, Paz, € Bolaños, 2010; Reyes et a/., 2011). 


Para los datos experimentales del maíz, el patrón lineal del 
espacio ao-bo de la Figura 2b representa la fase reproductiva, donde 
las propiedades ópticas de las hojas (mezcla hojas y órganos 
reproductivos) y/o las geométricas (distribución angular de las hojas y 
órganos reproductivos) muestran cambios. Así, el patrón temporal del 
TAF se muestra en la Figura 3b, donde después de la fase expo-lineal 
está una fase reproductiva antes de alcanzar el IAF máximo (IAFmx). 


Desde el punto de vista de las fases de los cultivos en la etapa 
vegetativa del crecimiento, las diferencias entre los patrones 
mostrados en la Figura 2 son fundamentales para entenderlos errores 
asociados con los modelos radiativos y para el diseño de índices de 


vegetación generalizados. 
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El problema fundamental del espacio del R-IRC es que induce un 
cambio de pendiente en la curvatura del patrón entre ao-bo, que 
muchos IV la interpretan como una “saturación” de las bandas. En 
realidad, no hay tal saturación, sólo es consecuencia de que los 
patrones lineales de algunos índices difieren más y más de los valores 
del segmento, aproximadamente lineal, asociado con ao-bo cuando el 
TAF es mayor que 2-3. Esta saturación es consecuencia del formato de 
razones usado en la mayoría de los IV con patrones lineales en el 


espacio paramétrico ao-bo. 


Índices espectrales de la vegetación 


Uno de los índices espectrales de la vegetación más utilizados es el 
NDVI, el cual está definido como (Rouse, Haas, Schell, Deering, € 
Hanlan, 1974): 


IRC —R 


add IRC +R (3) 


que implica una relación en el espacio ao-bo de ao = O y bo = f(NDVI) 
(Paz et al., 2015); esto es, supone que todas las líneas iso-IAF tienen 
origen en (0, 0) en el espacio R-IRC, por lo que en el espacio ao-bo de 
la Figura 2 representa una línea recta vertical sobre el eje bo (Paz et 


al., 2015). De la Figura 1 y Figura 2 resulta evidente que esta hipótesis 
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es incorrecta tanto para el caso de las simulaciones radiativas como 


para los datos experimentales. 


Para evitarlas hipótesis del NDVI, Paz et al. (2007) diseñaron un 
IV que toma en cuenta los cambios de la pendiente bo en función de 


do: 


aaa ; l=c+ da, (4) 


bo+1 bo 


NDVIcp = 


donde c y d son dos constantes empíricas asociadas con el patrón 
mostrado en la Figura 2b. 


La Figura 4 muestra los resultados de aplicar el NDVI y el NDVIgp 
al experimento de maíz (cuatro suelos de fondo con diferentes 
propiedades ópticas) presentado en la Figura 1b. Aparte de las 
variaciones de los valores del NDVI para el caso de suelos diferentes, 
este se satura alrededor de un JAF de dos, a diferencia del NDVI cp. 
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Figura 4. Resultados de aplicar índices espectrales al experimento 
del maíz: a) NDVI cp y b) NDVI. 


Tecnología y ciencias del agua, ISSN 2007-2422, 13(3), 374-424. DOI: 10.24850/j-tyca-2022-03-08 


2022, Instituto Mexicano de Tecnología del Agua 
Tecnología y Open Access bajo la licencia CC BY-NC-SA 4.0 
Ciencias¿Agua (https://creativecommons.org/licenses/by-nc- 
sa/4.0/) 


Paz et al. (2014 y 2015) han revisado la mayoría de los IV 
publicados en la literatura y concluyeron que ninguno aproxima en 
forma correcta el patrón expo-lineal definido en la Figura 2, incluido el 
NDVIcp, el cual fue diseñado considerando sólo el patrón exponencial 
y su transición en el espacio ao-bo. Para aproximar en forma completa 
el patrón expo-lineal entre ao y bo (Figura 2), Romero et a/. (2009) 
usaron una transformación de los espacios espectrales para el diseño 
del IV_CIMAS, el cual fue planteado para aproximar el patrón lineal 
hasta la transición al exponencial (lo contrario al NDVIcp). Bajo esta 
perspectiva, el IV_CIMAS está definido en un formato similaral NDVIcp 
(en función de las pendientes bo; Ecuación (2)), donde se tiene 
IV_CIMAS = NDVIcp cuando el valor de bo está debajo de un umbral 
dado e IV_CIMAS se obtiene del espacio espectral transformado 
cuando bo está arriba del umbral (Romero et a/., 2009; Paz et al., 
2015). 


Índices espectrales de la vegetación y Kc 


La estimación de la evapotranspiración actual del cultivo ha sido 
propuesta en función del coeficiente del cultivo K. promedio 
(Doorenbos 8 Pruitt, 1977) y de la evapotranspiración de referencia 
ETo: 


ETc = KcETo (5) 
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La relación (5) se refiere a condiciones promedio, por lo que 


Wright (1982) propuso el método dual de estimación de Kc: 


Kc = Kcb + Ke (6) 


donde K¿» es el coeficiente del cultivo de base (sin estrés) y Ke es el 


coeficiente asociado con la evaporación del suelo. 


En el caso de cultivos sujetos a estrés hídrico, el coeficiente del 
cultivo de base ajustado, Kcba, está dado por: 


Kcba = (Kcb)(Ks) (7) 


El método dual y su ajuste por estrés hídrico forman parte de la 
metodología actual propuesta por la FAO (Allen et a/., 1998), que usa 
un cultivo de referencia virtual, y es la base de las aproximaciones a la 


estimación de ET¿ usando el coeficiente del cultivo. 


Para entender las limitaciones de las relaciones entre los 
diferentes IV y el Kc (Kcb, Ks y Ke), la Figura 5 muestra el caso de las 
simulaciones radiativas del maíz para la condición de un suelo seco y 


uno húmedo del cultivo con y sin estrés hídrico (Reyes et a/., 2011). 
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Maíz - Simulación 


Sin estrés hídrico 


IRC (%) 


Línea del Suelo (IRC= as + bsR) 


Suelo húmedo 


0 5 10 15 20 
R (%) 


Figura 5. Simulaciones radiativas del cultivo de maíz bajo dos 
condiciones de humedad del suelo y de estrés hídrico. 


La humedad del suelo (K+) se manifiesta en curvas ¡so-suelo 
diferentes (R e IRC disminuyen sus valores) y el estrés hídrico (Ks) en 
curvas iso-IAF diferentes (menor pendiente para el caso de las líneas 
iso-TIAF). Ahora bien, considerando que la gran mayoría de los IV tratan 
de aproximar la pendiente (o la intersección) de las líneas iso-IAF (Paz 
et al., 2007; Paz et al., 2015), es claro de la Figura 5 que un IV 
“perfecto” (estimación directa, y sin error, de la pendiente) no puede 


diferenciar entre Kc» y Ke (o Kc), ya que la pendiente bo (o intersección 
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ao) de las líneas iso-IAF no contiene información del contenido de 
humedad del suelo. Así, el uso de un IV que aproxima bo, O ao, genera 
relaciones similares para Kcb, Kcóh + Ke O Kc. Esto impone restricciones 
en el uso de los IV para estimar el estado de la humedad del suelo, por 
lo que es necesario usar las líneas iso-suelo para este fin, a nivel de 
espacios transformados de las reflectancias (Paz et a/l., 2009) o con 


índices de vegetación (Paz, Reyes, € Medrano, 2011). 


Medios radiativos equivalentes 


En una parcela de cultivo, el IAF estimado experimentalmente 
representa uno global: 


IAF, - (2 foliar promedio ) pe de >) (8) 


planta m? de la parcela 


Sin considerar efectos de borde, el JAF; o local (plantas 
individuales) está relacionado con el global como: 


IAF, =1AF, fv (9) 


donde fv es la fracción de vegetación y está dado en proporción de O 
a 1. Así, la fracción de cobertura del suelo es fs = 1 - fv. 
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En la Figura 6 se muestra el IAF;, donde la fracción de huecos (fh) 
dentro del límite del follaje (proyección horizontal) estará dado por la 


Ley de Beer-Lambert para hojas no negras (Goudriaan € van Laar, 
1994): 


Área 1. huecos 
foliar 4 


Figura 6. Área foliar y huecos para el caso de un medio local o 
planta individual. 


f, =1—T,; T,= exp[(1 — w)?Kx14F,) (10) 
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donde T;¡es la transmitancia (local) de la radianza directa solar, w = p 
+ Tes el albedo foliar (suma de la reflectancia y transmitancia foliar, 
respectivamente) y K es un coeficiente de extinción total. Nótese que 
fv + fhn, ya que fv se refiere a la fracción de vegetación fuera de los 
límites de las plantas individuales, aunque en la práctica (fotografías 
digitales) fv incorpora a fh. 


En la Figura 7 se muestra el caso de un cultivo sembrado en 
surcos, donde al principio de la etapa de crecimiento se observa un 
patrón de plantas individuales, donde las plantas no se tocan entre sí 
(caso 3-D). Al avanzar el crecimiento, las plantas en un mismo surco 
se tocan entre sí y forman un cultivo en hileras o surcos (caso 2-D). 
Cuando las plantas de diferentes surcos se tocan, tenemos la situación 
de cobertura total del suelo (caso 1-D). Esta dinámica de crecimiento 
requiere de una estrategia genérica de modelación que tome en 
cuenta, razonablemente, la complejidad geométrica de los cultivos en 
parcelas y la vegetación natural. 


o0o0S 
O 


pea) 
e 


(a) (b) 


Figura 7. Geometría asociada con un cultivo sembrado en surcos: 
(a) suelo desnudo; (b) plantas individuales sin traslape; (c) plantas 
traslapadas en surcos, y (d) vegetación densa traslapada. 
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De la discusión anterior, un medio equivalente (términos físicos) 


es aquel que cumple (Paz et a/., 2013): 


_ IAF) 


IAF, == (11) 


Es decir, para una medición de JAFa hay un número infinito de 
combinaciones equivalentes de JAF, y fv que satisfacen la Ecuación 
(11). 


Aunque un medio radiativo equivalente contiene la misma 
cantidad de fitoelementos, las reflectancias pueden ser diferentes 
dependiendo de la cantidad de suelo expuesto y su condición de 
humedad. Así, los IV estimados y su relación con fv e IAFG es 
dependiente de la configuración geométrica del elemento suelo y 
vegetación en una parcela, así como de las propiedades ópticas 
asociadas. 


La discusión de medios radiativos equivalente, Ecuación (11), es 
importante para entender las limitaciones del uso de fv o TAFG para 
estimar en forma indirecta Kc (o Kcb), tal como se propone en FAO-56 
(Allen et a/., 1998). En este método se usa un valor de fv = 10 % para 
definir el inicio de la etapa vegetativa del crecimiento, cuyo final está 
dado por un fv de 70-80 % o un JAF de 2-3 (Allen eta/., 1998). De la 
Ecuación (11) es claro que es necesario definir un IAF; (fv e IAFG) para 
poder asociar en forma única la geometría y cantidad de vegetación 
(estomas) en una parcela, pues la definición de un solo valor, fv o IAFG, 


no es suficiente para caracterizarel medio (Odi-Lara, Paz-Pellat, López- 
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Urrea, € González-Piqueras, 2013). Este problema se manifiesta en la 
definición de la “cobertura efectiva” (fin de la etapa vegetativa), donde 
se le asocian valores de fv o IAF fijos (Allen et a/., 1998). Esto implica 
un cultivo estandarizado en relación con su geometría de siembra, y 
no capta la esencia de la arquitectura de las plantas y la distribución 
geométrica en las parcelas. Así, es posible encontrar que la cobertura 
efectiva se puede dar en valor de fv o JAFG, que pueden diferir 
sustancialmente de los recomendados por FAO-56 (Allen et a/., 1998). 


Materiales y métodos 


Se analizaron dos cultivos con densidad media de follaje, maíz (Zea 
mays L.) var. dracma 700 (Calera, González-Piqueras, 8. Meliá, 2004; 
González-Piqueras, Calera, Gilabert, Cuesta, 8 De-la- Cruz, 2004) y 
trigo (Triticum aestivum L.) var. estero (López-Urrea, Montoro, 
González-Piqueras, López-Fuster, 8 Fereres, 2009a). Se midieron las 
variables IAF, fv, Bm (biomasa) y h; además de medidas espectrales 


y evapotranspiración del cultivo en un lisímetro de pesada. 


Los experimentos se realizaron en el Instituto Técnico 
Agronómico Provincial de Albacete, España (39% 14' N, 29 5' O, 695 
msnm) en los años 2001 y 2003 para maíz y trigo, respectivamente. 
El clima es continental semiárido, con una precipitación anual promedio 


de 320 mm concentrada en primavera y otoño. Las temperaturas 
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media, máxima y mínima son 13.7, 24.0 y 4.5 %C, respectivamente 
(López-Urrea et al., 2006a; López-Urrea et a/., 2009a). El suelo fue 
clasificado como petrocalcic calcixerepts (Soil Survey Staff, 2006). La 
profundidad media es de 40 cm, limitado por un horizonte petrocálcico 
que se encuentra más o menos fragmentado. La textura es franco - 
arcillo-limosa, con un 13.42 % de arena, un 48.89 % de limo y un 
37.69 % de arcilla; el pH es básico. El suelo es pobre en materia 
orgánica y en nitrógeno total, con un alto contenido en caliza activa y 
potasio. El contenido de agua a capacidad de campo (CC) fue de 0.34 
y en punto de marchitez permanente (PMP) de 0.21 m3 m3, En el caso 
del trigo, la humedad aprovechable (HA) estimada para una zona de 
raíces de 35 cm fue de 45.5 mm y para la calendarización del riego se 
utilizó un nivel de agotamiento permisible de 0.55 recomendado por 
Allen et al. (1998) en trigo para valores de ET¿de 5 mm d*!, 


Los cultivos crecieron en condiciones óptimas, sin restricciones 
de agua y nutrientes. El trigo se fertilizó después de la siembra con 28 
kg hai de N, 112 kg ha"! de P205 y 70 kg ha? de K20; en la etapa 
vegetativa se aplicó fertilizante líquido (170 kg ha"? de N). El sistema 
de riego de la parcela es por aspersión mediante cobertura total 
enterrada. Se hizo lo posible para que el cultivo del lisímetro 
mantuviera la misma tasa de crecimiento y densidad de plantas que la 


parcela de protección. 


Para determinar la evapotranspiración del cultivo se utilizó un 
lisímetro de pesada con registro continuo de datos que se encuentra 
en el centro de una parcela de 100 x 100 m (López-Urrea et al., 
2006a). El lisímetro mide 2.7 m de largo, 2.3 m de ancho y 1.7 m de 


profundidad, con una masa total de 14.5 mg. La celda de carga tiene 
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una resolución de 0.25 kg, que equivale a 0.04 mm de agua. La medida 
de evapotranspiración se estimó a partir de los datos de peso del 


lisímetro, de la forma: 


ET=P+R-D+24w (12) 


donde P es precipitación, R es el riego, D el drenaje y 4w corresponde 


a las variaciones en el contenido de agua en la masa de suelo. 


Se trabajó con los datos horarios para obtener el dato diario de 
agua evaporada. Cuando no se pudo obtener el valor diario, por 
diferentes motivos (lluvia, riego, corte del césped), se hizo un 
promedio entre el valor anterior y siguiente a ese día. Otras fuentes de 
error que condujeron a la eliminación de datos fueron la verificación 
del peso, calibraciones por pérdida de contacto en la celda de carga, 
baja de voltaje en el sistema, pérdida de contacto del cable a tierra y 
falla en el controlador de dato (López-Urrea et a/., 2006a). 


La ETo se calculó mediante la ecuación FAO-56 de Penman- 
Monteith (Allen et a/., 1998), que ha tenido mejores resultados para la 
zona en estudios previos (López-Urrea et a/., 2006a; López-Urrea et 
al., 2006b). Los datos climáticos se obtuvieron de una estación 
meteorológica situada en la finca experimental. Los datos faltantes de 
radiación neta, velocidad del viento y humedad relativa se estimaron 
usando el dato previo. El paso de integración fue de una hora y a partir 
de estos datos se obtuvo el dato diario. 


Las mediciones espectrales se hicieron con un radiómetro GER 


3700 (Spectra Vista Corporation, NY, USA), que tiene un intervalo 
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espectral de 300 a 2500 nm, con una resolución espectral máxima de 
1.5 nm. La altura de observación de la cubierta en el caso de trigo fue 
1.80 m y en el maíz 5 m que permite un diámetro de observación de 


21 y 58 cm, respectivamente, para un campo de visión (FOV) de 6.60, 


Se hizo un seguimiento periódico de cada cultivo (entre 7 y 10 
días). En maíz se midió un transecto de seis puntos dentro y fuera del 
lisímetro para asegurar su representatividad respecto a la parcela, 
mientras que en trigo se hizo un barrido de cinco medidas sobre el 
lisímetro. Las medidas se hicieron en horario de máxima elevación 
solar y en días claros (sin nubes), manteniendo la observación a nadir 
(Milton, 1987). Se tomó la reflectancia de un panel de referencia 
Spectralon (LabsphereWR) antes de cada medida de la superficie para 
evitar que la configuración de iluminación variara (McCoy, 2005). 


Se registró la fenología de cada cultivo. En la Tabla 1 se detallan 
las fases de desarrollo. 


Tabla 1. Fenología de maíz y trigo. 


A Y AA 


Día Día 
juliano Fenología juliano Fenología 


pue e 
Mis ici ¿MUI 


Tres hojas 


Tecnología y ciencias del agua, ISSN 2007-2422, 13(3), 374-424. DOI: 10.24850/j-tyca-2022-03-08 


2022, Instituto Mexicano de Tecnología del Agua 
Tecnología y Open Access bajo la licencia CC BY-NC-SA 4.0 
Ciencias¿Agua (https://creativecommons.org/licenses/by-nc- 
sa/4.0/) 


Enderezamiento-un 
vii nudo 
Vvi14 


201 


Flor masculina Emergencia espiga 


5 
3 
Pe 


20 
Nic Sr Li 

219 Fecundación 148 Madurez lechosa 

225 o 000 608 Madurez pastosa 
A 

249 Pastoso 


clk O AS IS 


Maduración 


fisiológica 


La fracción de cobertura fv se estimó mediante la técnica de 
clasificación supervisada (máxima probabilidad) de fotografías digitales 
tomadas a nadir de la cubierta. La clasificación se hizo con el programa 
ENVIMR, distinguiendo entre vegetación verde y suelo. 


Para calcular el índice de área foliar, JAF, se recogió una muestra 
de 1 m2 de la parcela de protección y se midió una cara de las hojas 
con un equipo LI CORMR LAI 3000. 
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De acuerdo con las recomendaciones de FAO, se midió la altura 
del cultivo para adaptar los valores del coeficiente de cultivo Kc a lo 
largo del ciclo de crecimiento. La altura del cultivo se midió desde la 
base hasta el elemento más alto de la planta, seleccionando de forma 
aleatoria las plantas dentro de la parcela de protección. 


Resultados y discusión 


La estimación de Kc (Kc» + Ke) O Kecb (no distinguibles en términos de 
los IV) a nivel de predio o parcela agrícola requiere que el efecto del 
suelo (humedad) sea minimizado. En la perspectiva de un esquema de 
riego de precisión es importante conocer la variabilidad espacial de K. 


(ET¿) para hacer un uso eficiente del recurso agua. 


La Figura 8 muestra las estimaciones del NDVI para los diferentes 
días julianos, donde es evidente que este índice no reduce 
adecuadamente los efectos del suelo y el utilizar un promedio del NDVI 
es una decisión arbitraria, pues para un día particular, la condición de 
humedad del suelo en el tiempo de las mediciones permanece más o 
menos constante, por lo que las variaciones son atribuibles solamente 
a la incapacidad del NDVI de minimizar el efecto de suelo y/o 
variaciones en la condición del cultivo. Por otro lado, el NDVI muestra 
problemas de saturación (no varía en el tiempo) antes de que el cultivo 


alcance su JAF o fv máximo, limitando su uso sólo a la parte 
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exponencial de la etapa de crecimiento vegetativo, sin considerar la 
fase lineal. Los usos y costumbres de usar el NDVI para estimar Kc 
(Castañeda-Ibáñez et al., 2015; Palacios-Vélez, Palacios-Sánchez, 8 
Espinosa-Espinosa, 2018) no se justifican, ya que no es suficiente tener 
relaciones empíricas estadísticamente aceptables entre un IV y el Kc, 
es necesario tener las respuestas correctas por las razones correctas 
(Kirchner, 2006) para su uso. El no entender conceptos de 
transferencia radiativa asociados con los IV limita fuertemente los 
enfoques empíricos, justificándolos sólo en ajustes estadísticos 
arbitrarios. Debido al uso de la banda fotosintética (R) y la no 
fotosintética (IRC), se espera que cualquier IV aproxime, aunque sea 
de forma burda, el crecimiento de un cultivo y, por lo tanto, la 


evolución del Kc. 


1.0 
0.9 


: El 


0.5 
0.4 
0.3 
0.2 
0.1 
0.0 
100 120 140 160 180 200 220 240 260 280 300 


Día Juliano 


NDVI 


Figura 8. Variación del NDVI en el tiempo para las mediciones en el 
lisímetro para el maíz. 
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Las mayores vanaciones al inicio y final del ciclo del crecimiento 
del maíz, mostradas en la Figura 8, reflejan la condición de vegetación 
con baja cobertura del suelo, por lo que el efecto de este último 
predomina. 


Para evitar el problema de aproximar las pendientes de las líneas 
iso-IAF en el espacio del R-IRC, se estimaron directamente de los datos 
espectrales (sin efecto del suelo). Con las pendientes bo estimadas para 
cada día Juliano del ciclo de crecimiento del maíz y trigo, se estimó el 
IV_CIMAS, el cual representa en forma completa todo el ciclo. 
Asimismo, de las mediciones de K¿ en el lisímetro, se revisaron para 
seleccionar valores de Kc» (los valores más bajos) del patrón temporal 
de Kc. La Figura 9 muestra los patrones temporales de Kc«p e IV_CIMAS 
para el cultivo de maíz y trigo, donde se observa que hay un desfase 
en tiempo entre el crecimiento biofísico, representado porel IV_CIMAS, 
ya que hay una delación lineal entre este índice y el IAF (Paz et al., 
2007; Romero et al., 2009). 
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Figura 9. Patrones temporales de Kc» e IV_CIMAS: (a) maíz y (b) 
trigo. 


La Figura 10 ilustra las relaciones lineales entre el IV_CIMAS y 
Kcb, donde resulta claro que hay una relación para la etapa vegetativa 
y otra para la de senescencia, con valores constantes (etapa de 
meseta) de Kc» entre las dos etapas (cuadrados rellenos en la Figura 
10). Los parámetros de las relaciones lineales son diferentes entre la 


etapa vegetativa y la de senescencia, porlo que no se puede usar una 
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relación lineal única, como se ha propuesto por diferentes autores 
(Bausch, 1993; González-Piqueras et al., 2004; López-Urrea, De- 
Santa-Olalla, Fabeiro, Moratalla, 2009b; Palacios-Vélez et a/., 2018). 


2 
o 
sE 
0.00 0.20 0.40 0.60 0.80 1.00 1.20 
IV_CIMAS 
* Etapa Vegetativa 
o Etapa de Senescencia 
E] 
o 
pr 


y = 1.642x + 0.100 
R? = 0.996 


y = 1.7389x - 0.7711 
R? =0.9711 


0.00 0.20 0.40 0.60 0.80 1.00 1.20 


Figura 10. Relaciones entre el IV_CIMAS y el Kcb: (a) maíz y (b) 
trigo. 
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La Figura 11 ¡ilustra una relación lineal entre el IV_CIMAS y el IAF 
para el cultivo de maíz y muestra relaciones similares al caso de Kcb. 
Así, la parametrización de K¿ o Kepen función de fv o IAF es equivalente 
a la del IV_CIMAS. Por tanto, la parametrización de K¿« o Ke» como una 
función de fv o TAF es equivalente a las del IV_CIMAS. Dados los 
problemas de usar sólo IAF (o fv), las relaciones con Kc O Kcp varían de 
acuerdo con las condiciones locales; es decir, los parámetros de las 
relaciones lineales dependen de la geometría de siembra y desarrollo 
de los cultivos. 


6 
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IAF 
(99 


y = 19.267x - 5.6165 
R? = 0.9707 


y = 8.3312x - 3.7129 
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Figura 11. Relaciones entre el IV_CIMAS y fv e TAF. 


Del análisis de la problemática asociada con el uso de fv y IAF en 
forma separada, la definición del inicio y terminación (cobertura 
efectiva) de la etapa vegetativa es necesaria para parametrizar la 
relación con el IV_CIMAS. En este sentido se requiere una estimación 


directa del inicio y terminación usando los patrones temporales del 
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IV_CIMAS, que dada la relación lineal con el IAF, el modelo expo-lineal 
de la Ecuación (1) puede usarse para este fin (Odi-Lara et a/., 2013). 
Así, definiendo el inicio de la etapa vegetativa en el tiempo tr y el fin 
en tmax, podemos usar los valores del IV_CIMAS en esos tiempos y 
definir el Kc», O Kc, asociado con la etapa (datos de FAO-56 oO 
estimaciones con lisímetro), de tal manera que con el par de puntos 
poder determinar los parámetros de la relación lineal. Ahora bien, 
considerando que en términos del IV_CIMAS la etapa vegetativa está 
representada por la fase lineal (con un inicio de la transición 
exponencial a lineal) y que la relación temporal del K. o Kc» tiende a 
ser no lineal (exponencial) (Figura 9) podemos usar una 
transformación logarítmica-exponencial para hacer equivalentes 


ambos patrones (1 significa inicio y fes fin): 


= HDN. o =1n( Kcb, f) — fxIVemasf; Kcb = expl e + 


IVcimas ¡7!Vcimas y 


fxIVemas) (13) 


Usando la Ecuación (13), la Figura 12 muestra la relación entre 
los Kc» estimados y los medidos en el lisímetro para el maíz y el trigo, 
usando los valores medidos del lisímetro para parametrizar la Ecuación 
(13). 
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Figura 12. Relaciones el Kc» estimado y el medido: a) maíz y b) 
trigo. 


La relación entre los Ke» estimados, Ecuación (13), y los medidos 
mostrados en la Figura 12 muestra que la estimación con dos pares de 
datos para autocalibrar el IV_CIMAS da resultados buenos y puede 
usarse en términos operacionales. La relación entre el IV-CIMAS y Ku 
es similar para K¿, por lo que podemos usar los valores de Kc al inicio 


y término de la etapa vegetativa (fase lineal) en la Ecuación (13) y 
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obtener estimados de K¿ en la etapa. El valor de inicio (tiempo antes 
de tr) tiene un valor constante igual a Kcp,¡ O Kc,i. 


Algo similar al caso de la etapa vegetativa puede hacerse para la 
etapa de senescencia, definiendo así una metodología autocalibrable 
(cuadros de FAO-56 o mediciones de campo) para definir la relación 
entre el IV_CIMAS, y Kcb O Kc. 


La ventaja de la metodología propuesta es que es independiente 
de efectos atmosféricos, pues se calibra con las estimaciones locales 
del IV_CIMAS y los valores de Kc» O K¿ usados como inicio y terminación 
de una etapa. Bajo la consideración de efectos relativamente similares 
en un ciclo de crecimiento, o bajo la consideración de su minimización, 
las estimaciones resultan adecuadas para ser usadas en aplicaciones 
en tiempo cuasi-real, bajo las hipótesis o restricciones asociadas con 
la metodología FAO-56. 


Relación IV_CIMAS con la biomasa aérea 


La Figura 13 muestra la relación entre la biomasa aérea seca (Bm) y 
el IV_CIMAS para el maíz y trigo, donde para la etapa vegetativa (y la 
de la meseta) del IV_CIMAS se obtiene una relación exponencial con 
la Bm. Algo similar sucede para el caso de la etapa de senescencia. 
Esto es, usando IV_CIMAS versus In (Bm) se obtienen dos relaciones 


lineales que se interceptan en el inicio de la senescencia. 
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Figura 13. Relaciones entre el IV_CIMAS y la Bm seca: a) maíz y b) 
trigo. 
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Conclusiones 


Este trabajo describe los alcances y limitaciones del uso de Índices 
espectrales de la vegetación (IV) para la estimación de los coeficientes 
del cultivo (Kc O Kc») de la metodología FAO-56. 


Considerando que los IV tratan de aproximar la pendiente bo, o 
intersección ao, de las líneas iso-IAF de los cultivos, se aclaró que la 
relación entre los IV, y Kc» O Kc son equivalentes, pues por definición 
las líneas iso-TIAF no pueden estimar la humedad del suelo. Para esto 
último se requiere el uso de las líneas iso-suelo, tal como lo discuten 
Paz et al. (2009). 


Bajo la consideración de medios radiativos equivalentes 
(espacios espectrales), el uso de la fracción de cobertura (fv) o el índice 
de área foliar global (IAF) en forma independiente genera soluciones 
múltiples a la definición de los inicios y finales de las etapas del K¿ de 
acuerdo con FAO-56. Para parametrizar estos tiempos es necesania la 
consideración conjunta del IAF y fv para que refleje la geometría de 


siembra de los cultivos. 


Usando el formato genérico del índice IV_CIMAS, en conjunto con 
estimaciones directas de las pendientes de las líneas iso-IAF en dos 
experimentos (maíz y trigo), se obtuvieron relaciones lineales entre el 
Kcp y el IV_CIMAS para la etapa vegetativa y la de senescencia. Esto 
implica que no puede usarse una relación única entre los IV, y Kc O Kcb. 
De esta restricción, que requiere de calibraciones locales con 


mediciones de campo, se introdujo una metodología para la estimación 
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directa de la relación del IV_CIMAS, y Kcp O Kc para ser aplicada en 
términos operacionales. 


Los análisis mostrados en este trabajo definen las limitaciones en 
la estimación de K¿ O Kc» dada la limitación intrínseca del diseño 
convencional de los IV para aproximar las pendientes de las líneas iso- 
TAF. Así, bajo esta consideración, no es posible utilizar los IV para 
aproximar Ke (humedad del suelo). No obstante lo anterior, la 
estimación de Kc» afectado por estrés hídrico (evaluado de alguna 
manera), ¡.e., Kc»*Ks, puede realizarse en forma similar al Kc» sin 
estrés, siempre y cuando se tengan los valores iniciales y finales de 
Kce»*Ks asociados con las etapas de la metodología FAO-56. Este 


aspecto merece una mayor investigación. 
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